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Among the various techniques available, ultrasonic Lamb waves offer a convenient method of 
evaluating composite materials. Since the Lamb wave velocity depends on the elastic properties of 
a structure, an effective tool exists to monitor damage in composites by measuring the velocity of 
these waves. Lamb wave measurements can propagate over long distances and are sensitive to the 
desired in-plane elastic properties of the material. This paper describes two studies which monitor 
fatigue damage and two studies which monitor thermal damage in composites using Lamb waves. 
In the fatigue studies, the Lamb wave velocity is compared to modulus measurements obtained 
using strain gage measurements in the first experiment and the velocity is monitored along with the 
crack density in the second. In the thermal damage studies, one examines samples which were 
exposed to varying temperatures for a three minute duration and the second includes rapid thermal 
damage in composites by intense laser beams. In all studies, the Lamb wave velocity is 
demonstrated to be an excellent method to monitor damage in composites. © 1998 Acoustical 
Society of America. [S0001-4966(98)02705-2] 

PACS numbers: 43.35.Zc, 43.40.Le, 43.20.Mv [HEB] 


INTRODUCTION 

Composite materials are being more widely used today 
by aerospace, automotive, sports equipment, and a number of 
other commercial industries because of their advantages over 
conventional metals. Composites have a high strength-to- 
weight ratio and can be constructed to meet specific design 
needs. In addition to satisfying the design specifications, 
composites offer weight savings over traditional metals with- 
out sacrificing strength. With the advantages of composite 
materials, their use is continually increasing and becoming 
more widespread to a variety of fields. 1-3 

O’Brien et al. 4 and Charewicz and Daniel 5 have exam- 
ined the fatigue behavior of laminates and proposed lifetime 
prediction models based on fatigue factors such as increased 
matrix cracking, delamination growth, and stiffness loss. 
Ogin et al. 6 ' 1 and Lim and Hong 8 have modeled the stiffness 
reduction in composites due to the effects of matrix cracking. 
Laws and Dvorak 9 have investigated the progression of ma- 
trix cracking as well as stiffness loss due to transverse crack- 
ing. 

In addition to fatigue damage, thermal degradation of 
composites is of concern. Street et al . 10 have examined how 
composites respond to thermal damage of temperatures up to 
350 °C for short durations, as well as exposure times as long 
as 30 min for temperatures of 225 °C. Herakovich and 
Hyer" and Ermer and Mall 12 have studied the effects on 
composites which were thermally cycled. The temperature 
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range used was - 157 °C-+121 °C in Ref. 11 and 149 °C- 
427 °C in Ref. 12. Stansfield and Pritchard 13 have subjected 
composites to rapid changes in temperature. These fluctua- 
tions, similar to the flight profile of supersonic aircraft, can 
lead to moisture absorption and matrix cracking. 

Traditional nondestructive evaluation (NDE) techniques 
used on metals are often not easily adaptable to composites 
due to the inhomogeneous and anisotropic nature of these 
materials. Finding an effective means of nondestructively 
monitoring damage in composites is a challenging task, but it 
is extremely important to the safety and reliability of such 
structures. Several techniques are available to nondestruc- 
tively evaluate damage in materials. 14-16 Salkind 14 and 
Stinchcomb 15 detail various methods used for damage detec- 
tion in composites such as ultrasonics, x-ray radiography, 
thermography, and holography. In addition, Daniel et al. 16 
examined the use of ultrasonic backscatter to monitor crack 
density and compared the results to those obtained using 
x-ray radiography. 

Among the various techniques available, ultrasonic 
Lamb waves offer a convenient method of evaluating these 
composite materials. As a material is damaged, the elastic 
parameters of the structure change. Since the Lamb wave 
velocity depends on these properties, an effective tool exists 
to monitor damage in composites by measuring the velocity 
of these waves. Additionally, Lamb wave measurements are 
better than conventional through-the-thickness ultrasonic 
measurements because they can propagate over long dis- 
tances and are sensitive to the desired in-plane elastic prop- 
erties of the material. 



Studies have been conducted which show a reduction in 
Lamb wave velocity due to a loss of stiffness caused by 
matrix cracking. 17-19 All of the authors noted that the in- 
plane properties of a composite are affected by cracking and 
that they can be monitored through the use of Lamb waves, 
which travel perpendicular to the crack direction. Dayal and 
Kinra 17 showed that for through-the-thickness measure- 
ments, the attenuation is sensitive to cracking, while the ve- 
locity is not. For Lamb wave propagation in the plane of the 
plate, both wave speed and attenuation were found to be 
sensitive to cracking. Tang and Henneke 18 noted that Lamb 
waves provide information about the in-plane properties of a 
plate, and this type of measurement is more useful due to the 
fact that composites are commonly designed to carry in- 
plane loads. Similarly, Dayal et al . 1 y noted that the Lamb 
wave interaction with cracks is much stronger in the plane of 
the plate and, thus, provides an effective method to detect 
damage due to transverse matrix cracking. 

Given the ability to nondestructively monitor materials 
using guided waves, studies have been conducted which use 
Lamb waves to monitor a variety of material properties in 
composite structures. Karim et al . 20 and Mai et al. 21 have 
used inversion techniques to determine the material param- 
eters of composites from experimental Lamb wave data. Ha- 
beger et al. 22 studied the use of plate waves to monitor the 
strength of paper for the purpose of monitoring quality con- 
trol on line for the paper industry. Several authors 23 25 have 
used Lamb waves to inspect adhesively bonded joints. Lamb 
wave techniques have also been used to study 
delaminations, porosity, and fiber misalignment. 
Recently, Lamb waves have even been used by Lee et a/. 31 
to measure the temperature of silicon wafers during process- 
ing and by Pei et al. 22 for in situ measurements of thin film 
thickness. Both of these techniques would be of interest to 
the semiconductor processing industry. 

While numerous studies exist which examine fatigue 
damage in composites using nondestructive techniques, few, 
if any, have been conducted which monitor thermal damage 
using ultrasonic nondestructive evaluation techniques. In the 
studies mentioned earlier, 10 13 the effects of thermal damage 
were measured using mainly mechanical tests and micro- 
scopic observation. This work will focused on the use of 
Lamb waves to nondestructively monitor fatigue as well as 
thermal damage in composite structures. 

The following sections describe two experimental stud- 
ies which monitor fatigue damage in composites using Lamb 
waves. In the first study, the Lamb wave velocity was com- 
pared to modulus measurements obtained using strain gage 
measurements. Additionally, experimental dispersion curves 
are compared for damaged and undamaged specimens. In the 
second study, the Lamb wave velocity was monitored along 
with the crack density. For both studies, the Lamb wave 
velocity measurements showed a significant decrease as the 
samples underwent increasing amounts of damage. 

In addition to the fatigue studies, two studies were also 
conducted which monitor thermal damage in composites us- 
ing the Lamb wave technique. The first study involves 
samples which were heated to varying temperatures for a 
three minute duration. The second study examined rapid 


thermal damage in composites. These samples were sub- 
jected to a very short duration (0.5 and 1.0 s) laser source. In 
both studies, the Lamb wave velocity was demonstrated to 
be an excellent method to monitor thermal damage in com- 
posite materials. 

I. FATIGUE DAMAGE 1 

In this experiment, studies have been conducted which 
monitor fatigue damage in composite samples using strain 
gage measurements as well as Lamb wave velocity measure- 
ments. A description of the test samples is followed by the 
results of two different measurements of Lamb wave veloc- 
ity. The first technique is a contact measurement done at a 
single frequency, while the second involved an immersion 
study of Lamb waves in which dispersion curves were ob- 
tained. The results of the Lamb wave method for monitoring 
fatigue damage are compared to the damage progression 
measured using strain gages. 

The composite samples studied were AS4/3501-6 
graphite/epoxy with a stacking sequence of [0/90 3 ] s . Two 
30.5-by 38.1 -cm plates were manufactured and C-scanned 
prior to being cut into specimens to check for any abnormali- 
ties. The scans revealed a 5%-7% level of porosity in the 
samples. The plates were then cut into 28-by 3.8-cm cou- 
pons. The finished coupons had an average thickness of 0.12 
cm. Two 0.635-cm strain gages were attached to each 
sample: one axial and one transverse. From the values mea- 
sured for the stress, axial strain, and transverse strain, one 
can obtain Young’s modulus and Poisson’s ratio for the 
sample. 

Before the fatiguing process, two samples were loaded 
quasi-statically to failure to obtain Young’s modulus, Pois- 
son’s ratio, and ultimate strength. Each of these coupons was 
loaded from 0 to 3250 microstrain and back twice. On the 
third run, the samples were loaded to failure. The modulus 
and Poisson’s ratio were obtained by averaging the values 
measured for each of the three load cycles described above. 
The ultimate strength was defined as the load at which the 
samples failed. 

The specimens were subjected to tension-tension fa- 
tigue in a 55-kip capacity load frame. The samples were 
fatigued at a frequency of 10 Hz and at an R value (mini- 
mum load/maximum load) of 0.3. The upper load was taken 
to be 33% of the ultimate strength, which had a value of 160 
MPa. The specimens were successively raised to higher fa- 
tigue cycle values. The velocity of the lowest order symmet- 
ric Lamb mode was measured and the modulus obtained 
from strain gage measurements both before and after each 
cyclic loading. The samples were removed from the load 
frame at intermediate values of fatigue cycles in order to 
make the contact measurements. 

Lamb modes can be generated by placing a longitudinal 
contact transducer on a plate and driving it with a sine wave 
tone burst. Ditri et al 2 ' showed that periodic loading normal 
to the surface of a plate, as produced by a longitudinal trans- 
ducer, will excite Lamb waves at distances away from the 
source which are large compared to the wavelength of the 
Lamb mode. Lamb waves are generated in most materials at 
a distance of roughly 1 cm or less. Since measurements are 




FIG. 1 . Experimental arrangement for the contact measurement. 


normally done anywhere from 5 to 10 cm away from the 
source, this does not create any difficulties in generating and 
measuring Lamb modes using this technique. 

A Lamb wave was generated by using a 2-cycle, 0.5- 
MHz sine wave to excite the sending transducer. The signal 
was received at various positions, digitized at 10 MHz, and 
the waveforms saved on a computer. Both the sending and 
receiving transducers were broadband longitudinal wave 
transducers with a nominal center frequency of 1 .0 MHz and 
a diameter of 1 .27 cm. A coupling gel was used between the 
transducers and the plate. The experimental arrangement is 
shown in Fig. 1. 

The time differences were measured by imposing a com- 
putational delay on the first signal to overlap the signals 
received at greater distances. The small thickness of the 
plates combined with the low frequency used to generate the 
Lamb wave yield a frequency-thickness product of 
0.6MHz-mm. In this region, only the S 0 and A 0 modes 
propagate. Therefore, the leading part of the wave was iden- 
tified as the lowest order symmetric Lamb wave, which is 
not very dispersive. The trailing portion of the signals con- 
tained the dispersive antisymmetric wave as well as edge 
reflections. The distances were measured to an accuracy of 
0.01 cm using a sonic ranger. A least-squares fit from a plot 
of time and distance was performed to obtain the velocity of 
the prompt S Q mode. 

As the samples underwent increasing amounts of fa- 
tigue, an increasing number of cracks developed in the 90° 
plies. These cracks degrade the modulus. Photomicrographs 
at different magnifications were taken of the edge of a fa- 
tigued composite and are shown in Fig. 2. For the lowest 
magnification, the fibers from the 0° laminas are the vertical 
bright streaks on the outer edges and the six 90° laminas are 
what appear as dots in between the outer layers. The ends of 
the fibers for these laminas is seen quite clearly at the higher 
magnification shown in Fig. 2(d). The size of the crack 
shown is approximately 1 /rm. 

The composite samples exhibited a decrease in the ve- 
locity as the modulus degraded due to increasing fatigue. 
Both the modulus and velocity squared follow the same gen- 
eral trend with increasing fatigue cycles. This is shown sche- 
matically in Fig. 3. In the plot, the velocity was normalized 
to the velocity measured before fatiguing. The normalized 



FIG. 2. Photomicrographs of the edge of a composite sample: (a) without a 
crack and ( b ) — (d ) with a transverse matrix crack. The magnification for each 
was: (a) 5.5; (b) 5.5; (c) 1 1.0; and (d) 28.0. The thickness of the sample is 
0.12 cm. 

modulus was obtained in a similar manner. The velocity 
squared error bars are an estimate of the total error due to the 
uncertainties in the distance measurement and the digital sig- 
nal overlap used to obtain the time differences. Some of 
values for the three measurements on the individual samples 
varied as much as 8% -10% from the average. Using a simi- 
lar technique, Tang and Henneke 18 have reported deviations 
as great as 10% from the average of three measurements as 
well. 

The velocity of the lowest order symmetric mode de- 
pends on the square root of the modulus for isotropic plates 
and the square root of the in-plane stiffness values for uni- 
directional composites. 34 Since the modulus shown in the 
plot was obtained from strain gage measurements of the in- 
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Fatigue Cycles 

FIG. 3. Plot of normalized modulus (solid circles) and normalized velocity 
squared (open triangles) versus fatigue cycles for a fatigued sample. 
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FIG. 4. Plot of the lower modes of the experimental dispersion curves for a 
sample which is undamaged (solid circles) and one which has been fatigued 
10 6 cycles (open triangles). 

plane stiffness, the square of the velocity was chosen to plot 
against the number of fatigue cycles for convenience. 

In addition to the contact measurements, which were 
done at a single frequency, an immersion method was used 
to obtain dispersion curves for two samples. The experimen- 
tal arrangement for this type of measurement is detailed by 
Balasubramaniam et al. 35 In this method, the sample was 
immersed in water, which serves as the couplant between the 
transducer and the plate, and the same transducer was used to 
generate and receive the ultrasonic signal. An aluminum 
block served as a reflector for the input pulse as the specimen 
was rotated about the corner. When a Lamb wave is excited 
in the plate, part of the total energy is converted to a propa- 
gating Lamb mode. Due to this loss of energy, a minimum in 
the frequency spectrum of the reflected signal corresponds to 
the frequency at which a Lamb mode is propagating. 

The reflected waveform at each angle was saved for later 
analysis of the frequency content. The velocity can be deter- 
mined from Snell’s law using the fact that only a longitudi- 
nal wave is supported in water and that the Lamb mode is 
propagating in the plate. After determining the frequency of 
the minimum for each angle, the velocity verses the 
frequency-thickness product can be plotted to obtain the dis- 
persion curve. The frequency range is determined by the 
bandwidth and frequency of the transducer used. 

In this experiment, an immersion transducer with a 
nominal center frequency of 1.0 MHz and diameter of 2.86 
cm was used in a pulse-echo arrangement. The specimen was 
rotated 0.4° -60° in steps of 0.4° and the reflected waveform 
at each angle was saved for later analysis of the frequency 
content. After determining the frequency of the minimum for 
each angle, the velocity (from Snell’s law) verses the 
frequency-thickness product was plotted to obtain the dis- 
persion curves. This technique was used to obtain dispersion 
curves for an undamaged sample as well as one which had 
one million fatigue cycles. The results are shown in Fig. 4. 

The dispersion curves for the fatigued sample is clearly 
shifted from the undamaged sample for both of the lower 
modes. The curves have roughly a 7% difference in velocity 
for the symmetric mode. This is in good agreement with the 
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FIG. 5. Comparison of lower modes of numerical curves using full param- 
eters (solid circles) and with c 33 reduced by 5% (open triangles). 

contact result for the same sample, which had a 7.2% relative 
decrease in velocity after one million fatigue cycles. It is also 
noted in the figure that the A 0 mode velocity variation was 
much less than that of the S 0 mode. This can be understood 
by examining how the velocity relates to the modulus. The 
velocity changes as the square root of the stiffness for the 
symmetric mode, but as the fourth root of the stiffness for 
the antisymmetric mode. 34 Therefore, the A 0 mode velocity 
will change less than the S 0 mode velocity for a given de- 
crease in modulus. 

For propagation across the fibers, where the 1-axis is in 
the fiber direction and the 3-axis is perpendicular to the fi- 
bers, a decrease in the elastic stiffness parameter c 33 alters 
the dispersion curves. This is due to the fact that c 33 is the 
modulus across the fibers. The transverse matrix cracks (due 
to fatigue damage) should, therefore, lead to a decrease in the 
elastic stiffness constant c 33 . To examine the effect of the 
matrix cracking, a through-the-thickness finite element 
model, detailed in Dong and Huang 36 and Datta et al., 31 is 
used to generate dispersion curves in laminated composites. 
Numerical dispersion curves were calculated in which the 
value of c 33 was reduced by 5% in only these layers. A 
comparison was done for the numerical dispersion curves for 
the lower modes. The curves for the lower modes calculated 
using the full elastic parameters and those calculated using 
c 33 reduced by 5% in the 90° layers are shown in Fig. 5. 
When compared in the region where the experimental data 
was obtained (O.lS^fd^ 1.05 MHz -mm), the numerical re- 
sults (see Fig. 5) are in good agreement with the experimen- 
tal dispersion curves shown in Fig. 4. For the theoretical 
curves in this region, the S Q velocity shows roughly a 7% 
decrease and the A 0 mode variation is much less than the S 0 
velocity change. These are exactly the characteristics which 
were observed in the experimental curves. 

It has been noted by Charewicz and Daniel 5 that strain 
gages only measure a local change in modulus and do not 
measure a global change in the specimen. Since the strain 
gages used in this study had a gauge length of 0.635 cm, they 
may not give an accurate account of damage throughout the 
entire plate. This would be especially true at the lower fa- 
tigue levels where very few cracks were observed in the 
samples. As the samples were damaged further, the cracking 




became more uniform. Thus a Lamb wave technique at low 
frequencies used to track the modulus as a function of fa- 
tigue history may provide improved global interrogation over 
the propagation length of the wave. 

II. FATIGUE DAMAGE 2 

Additional composite samples of the same material, 
AS4/3501-6, and same architecture, [0/90^] s , were obtained 
in order to conduct a second fatigue damage study. The new 
composites were loaded to greater levels in order to induce 
damage more quickly in the samples. Two composite plates 
identical to those used in the previous study were obtained. 
The new composite plates were C-scanned prior to cutting 
them into coupons to check for porosity. The images were 
very uniform and, therefore, no large scale areas of porosity 
were expected. The plates were then cut into 18 3.8-by 28- 
cm coupons. The edges were polished and photomicrographs 
were taken. No appreciable amounts of porosity were found. 

Similar to the first study, the specimens have been sub- 
jected to tension-tension fatigue. The samples were fatigued 
at a frequency of 10 Hz and at an R value of 0.35. For this 
study, however, the upper load was taken to be 45% of the 
ultimate strength, which has a value of 200 MPa. The 
samples were fatigued for a range of 10" cycles (n = 3, 4, and 
5). As before, the samples were removed from the load frame 
at intermediate values of fatigue cycles and the velocity of 
the lowest order symmetric Lamb mode was measured. In 
addition to loading these samples at higher levels, the crack 
density was also monitored during the process. 

A Lamb wave was generated by using a 6-cycle, 0.5- 
MHz sine wave to excite a longitudinal wave transducer with 
a nominal center frequency of 0.5 MHz and a diameter of 
1.27 cm. The signal was received at various positions by an 
array of pinducers with active elements of 0.135 cm. The 
pinducers were placed in a holder at a fixed separation dis- 
tance of 1 .0 cm between each transducer. This receiving ar- 
rangement greatly reduces the error in the measurement of 
the velocity due to the uncertainty in the distance between 
the transducers. Since a differential measurement is being 
made, fixing the transducers in a holder virtually eliminates 
the distance error. A coupling gel was used between the 
transducers and the plate. 

The signals were digitized at 10 MHz, averaged, and the 
waveforms saved to the computer. In an effort to reduce the 
error incurred in the time measurement by human judgement 
of the overlap region, the time differences between the pin- 
ducers were measured by performing a cross correlation. The 
leading part of both waves, which is identified as the lowest 
order symmetric Lamb wave, was used in the correlation. 
The distances between the five pinducers in the array were 
held fixed at 1, 2, 3, and 4 cm, respectively. The velocity of 
the prompt S 0 mode was obtained by a least-squares fit from 
a plot of time and distance. 

To measure the crack density, a 2.5-cm section in the 
middle of the plate was marked. This was done in order to 
count the number of cracks occurring in the same 2.5-cm 
length in the sample. Keeping the same measurement area is 
especially important at the lower cycle levels where the 
cracking is not very uniform. The number of cracks was 


counted by examining the edges of the coupons under a mi- 
croscope. The velocity measurements described above were 
centered over the same region where the crack density was 
monitored. 

The influence of fatigue cycles on transverse matrix 
cracking and Lamb wave velocity was examined for five 
different samples. The velocity and crack density were ob- 
tained after each cyclic loading. Each velocity was obtained 
from the average of three measurements. The standard devia- 
tion for all of the samples was less than 5% from the average 
of the three measurements and most were within 2.5% of the 
average for each individual sample. Thus the reproducibility 
of the measurements is quite good and provides an improve- 
ment over the technique used in the first fatigue study. Also, 
these results are similar to those obtained by Tang and 
Henneke, 18 who report a variation of 1.5% using a similar 
technique. 

As expected, the crack density for these samples in- 
creased with increasing cycles and, as with the previous 
study, the velocity squared decreased with fatigue cycles. 
Since strain gages were not used on these samples, the de- 
crease in modulus will be estimated from the crack density. 
An expression which relates reduced modulus to crack den- 
sity is given in Caslini et al. 3S For a [0„,/90„] s laminate, the 
stiffness loss is given by 38 

E nE 2 tanh(A./2D) _1 

~E~ 0 ~ l + ~^E[ {X/2 D) ’ ^ 

where 

2 _ 3Gi 2 E 0 (n+m) 
h 2 E l E 2 n 2 m 

In the above equations, E and E 0 are the values of the modu- 
lus for the damaged and undamaged laminate, respectively. 
E | , E 2 , and G 12 are the moduli (longitudinal, transverse, 
and shear, respectively) for each individual unidirectional 
lamina in the plate. D is the crack density, h is the lamina 
thickness, and m and n are the number of 0° and 90° layers, 
respectively, in the half-thickness of the plate. 

For the case presented here, m= 1, n = 3, and h 
= 0.15 mm. The value for the undamaged modulus, E 0 , was 
obtained from the average of the four modulus measure- 
ments for the samples which were loaded to failure. The 
engineering constants for AS4/3501-6, are taken to be 1 

£j = 142.0, E 2 = 10.34, G 12 =7.17, 

where all values are listed in GPa. 

The reduced stiffness as a function of crack density was 
calculated from the parameters given above using Eqs. (1) 
and (2). Figure 6 shows a comparison of the calculated re- 
duced stiffness and the normalized velocity squared for all of 
the samples in this study. The velocity was normalized in a 
similar manner as described previously. Since the uncer- 
tainty in the distance and time measurements in this study is 
less than that of the earlier study, the velocity error bars are 
not as large as were previously shown. Due to the possibility 
of overlooking a small crack under the microscope, the error 
in calculation of the crack density is assumed to be a ±1 
error in the total number of cracks counted in the 2.5-cm 
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FIG. 6. Calculated normalized modulus (solid line) versus normalized ve- 
locity squared (solid circles) as a function of crack density for all samples. 

length. However, for clarity, the error bars for the crack den- 
sity have been omitted in the figure. Figure 6 shows a larger 
decrease for the normalized velocity squared than for the 
normalized modulus which was calculated from the crack 
density. Although not as pronounced as the results shown in 
Fig. 6, this behavior was evident in the previous study when 
the modulus was measured using strain gages (see Fig. 3). 

These results are also in agreement with those found by 
Dayal et a/., 19 who examined the modulus degradation of an 
AS4/3502, [0/903] s composite using velocity measurement 
of the lowest order symmetric Lamb mode. The architecture 
and the fibers (AS4) of the samples used by Dayal et al . 19 are 
identical to the ones studied here. The epoxy matrix used 
(3502) is slightly different than the one used here, but the 
material properties of the two samples are very close. 

Similar to the measurements done here, the tests by 
Dayal et al . 19 were also done at a frequency of 0.5 MHz. The 
velocity of the S 0 mode, measured using a different tech- 
nique, and the crack density, measured using edge replica- 
tion, were monitored as the samples underwent increasing 
damage. The modulus was calculated directly from the ve- 
locity measurement using the relation 19 

E = pv\ (3) 

where E is the modulus, p is the density, and v is the velocity 
of the S 0 Lamb mode. Given this relation, the normalized 
modulus will also be equal to the normalized velocity 
squared. The results from Dayal et al . 19 are presented, along 
with the results for a representative sample from this study, 
in Fig. 7. 

In the figure, the solid line represents the theoretical re- 
duced modulus previously calculated using the material pa- 
rameters for AS4/3501-6. However, as mention earlier, the 
parameters for this type of composite are extremely close to 
those of AS4/3502. Therefore, the theoretical curve shown in 
the figure is representative of the composites used in both 
studies. Another interesting feature in the figure involves the 
comparison of the two velocity measurements. Both sets of 
data show a larger decrease in normalized velocity squared 
than what is theoretically predicted for the modulus reduc- 
tion. Dayal et al . 19 reported similar results when the modulus 
reduction measured using the Lamb wave technique was 
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FIG. 7. A plot of the normalized velocity squared for a fatigued sample 
(solid circles) versus data from Dayal et al. (Ref. 19) (open triangles). Curve 
fits are shown for the data in this study (short dashes) and for the data of 
Dayal et al. (Ref. 19) (long dashes). The solid line is the theoretical modu- 
lus reduction described in the text. 

compared with the results of analytical modeling and static 
stiffness measurements conducted by other authors. Thus the 
Lamb wave velocity technique may provide a more sensitive 
method of damage detection than modulus reduction mea- 
surements. 

Presented in the second fatigue study are the results of a 
different technique than that which was used in the previous 
study for determining the Lamb wave velocity in composites. 
Due to a reduction in the uncertainty in the distance and time 
measurements, the error in the velocity measurements were 
found to be less using this technique than those found in the 
first study. The crack density was also monitored as a func- 
tion of fatigue cycles. A decrease in normalized velocity 
squared with increasing fatigue cycles was noted. This was 
expected and agrees with the results presented in the initial 
study. 

The main difference between this study and the initial 
fatigue study was the monitoring of the crack density with 
fatigue cycles. The crack density was observed to increase 
with an increasing number of cycles, as is to be expected, 
and the Lamb wave velocity exhibited a decrease with in- 
creasing crack density. This drop in velocity as a function of 
crack density was compared with the results of Dayal 
et al ,, 19 who saw a similar velocity decrease. 

As was mentioned before, strain gages only measure a 
local change in modulus and do not measure a global change 
in the specimen. The strain gages used in the previous study 
had a gauge length of 0.635 cm and therefore, may not have 
given an accurate account of damage throughout the entire 
plate. It was shown in Caslini et al , 38 that the crack density 
can be theoretically related to stiffness reduction. Since the 
crack density was monitored over a 2.5-cm length, this be- 
comes the effective gauge length for measuring the modulus 
reduction via the crack density. 

The modulus reduction, calculated from the crack den- 
sity, was compared to the results obtained using the Lamb 
wave velocity measurements. A more dramatic decrease was 
observed in the velocity measurements than in the theoretical 
modulus reduction. Thus a Lamb wave technique used to 
track damage in composites may provide improved interro- 
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FIG. 8. Plot of S 0 velocity for a range of samples exposed to different 
temperatures for durations of 3 min (solid circles) and 10 min (open tri- 
angles). 

gation of the effective stiffness of a structure over the propa- 
gation length of the wave. 

III. THERMAL DAMAGE 

Determining the response of composites to thermal dam- 
age is necessary for a complete understanding of the total use 
environment of these materials. Heat damaged samples, 
which simulate typical thermal damage on carbon/epoxy air- 
craft structures, were provided by McDonnell Douglas Aero- 
space. These 15.24-by 15.24-cm specimens were fabricated 
as 16-ply, uniaxial plates of AS4/977-3 and had an average 
thickness of 3.07 mm. All samples had received different 
levels of heat damage, had been C-scanned, and had Barcol 
hardness readings performed. The exposure temperatures 
ranged from 288 °C to 538 °C for a duration of 3 min and 
included a sample which was exposed to 316 °C for 10 min. 

Lamb waves were generated and received with the same 
techniques described in the second fatigue study. The only 
differences were that the transmitting transducer had a nomi- 
nal center frequency of 1 .0 MHz and the receiving transduc- 
ers had fixed separations of 4, 5, 6, and 7 cm. The signal was 
received at various positions and digitized at 10 MHz. The 
time differences were measured by imposing a computa- 
tional delay on the first signal to overlap the signals received 
at greater distances. A least-squares fit from a plot of time 
and distance was performed to obtain the velocity of the 
prompt S 0 mode. 

As in the fatigue studies, the velocity of the lowest order 
symmetric Lamb mode was measured. The results of the 
velocity measurements in the fiber direction for ten different 
thermally damaged samples are plotted in Fig. 8. The error 
bars are an estimate in the velocity error and are mainly due 
to uncertainty in the digital channel alignment. The decrease 
in velocity with increasing thermal damage is quite dramatic. 
The variation in velocity is on the order of 18% for the two 
extremes of temperature exposure. The correlation with ex- 
tended temperatures is excellent and provides a method to 
assess thermal damage. Although the velocity showed a large 
change with thermal damage, the Barcol hardness numbers 
provided by McDonnell Douglas showed little variation over 



posite. 

this temperature range. Therefore, measurement of the Lamb 
wave velocity could provide a better technique to assess ther- 
mal damage than conventional methods. 

A model calculation 36 ’ 37 to determine the dispersion 
curves for this material system was computed using engi- 
neering constants supplied by McDonnell Douglas. The four 
standard engineering values (E x > e 2 > g 12 , and v l2 ) can be 
used to calculate the corresponding five elastic stiffness pa- 
rameters (cn, c 13 , c 23 , c 33 , and c 55 ) using various 
techniques. 1 Since a value for G 23 was not available, the 
value of c 44 [where c 44 =(c 22 — c 23 )/2] must be estimated 
using common values for polymer matrix/graphite fiber 
composites. 39 The value of C44 typically ranges from 3.5 to 
4.0 GPa, so a median value of 3.75 GPa was chosen for the 
calculation. The density of the material was taken to be 
1560 kg/m 3 . A five layer model calculation of the dispersion 
curves for propagation in the fiber direction is shown in Fig. 
9. 

From the dispersion curve, one can see that the S Q mode 
is almost completely nondispersive for the lower frequency- 
thickness products. Also, the So velocity in this region is 
around 9 km/s, which is exactly the velocity that was mea- 
sured for the samples with little or no damage. Thus the 
model compares well with the experimentally determined ve- 
locity. 

Theoretical dispersion curves were also calculated for 
AS4/977-3 using reduced stiffness constants. The effect of 
reducing different c ,/s by 20% on the low frequency, non- 
dispersive region of the lowest order symmetric mode (0.0 
S/rfSl.OMHz'mm) was investigated. The only parameter 
which altered the velocity of the S 0 mode in this region for 
propagation along the fibers is c x x . 

Figure 10 shows a plot of the experimentally determined 
normalized velocity versus exposure temperature, where the 
normalized velocity is the velocity in the damaged sample 
divided by the velocity in an undamaged control sample. 
Also plotted is the normalized velocity of the S 0 mode with 
reduced values of C[[. Here, the normalized velocity is the 
velocity with c u reduced divided by the velocity with the 
full values of the elastic stiffness constants. Both velocities 
were taken at a constant frequency-thickness product of 
0.5 MHz -mm on the dispersion curve. Theoretical calcula- 
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FIG. 10. Plot of the normalized S 0 velocity of samples exposed to different 
temperatures. The duration times are 3 min (solid circles) and 10 min (open 
triangle). Also shown (dashed lines) is the theoretically determined normal- 
ized S 0 velocity with reduced values of c n , as described in the text. 

tions with reduced values of Cp show that a 25%-30% re- 
duction in this parameter lowered velocities of the S 0 mode 
by 15%. This value corresponds to the velocity reduction of 
the sample which saw an exposure temperature of 538 °C. 

IV. LASER DAMAGE 

The three laser damaged samples examined were fabri- 
cated as quasi-isotropic plates with a stacking sequence of 
[0/+45/— 45/90] s . The fiber/matrix used on two of the 
samples was IM7/SP500 and the remaining sample was IM7/ 
PEEK. All three samples were irradiated at 100 W/cm 2 . The 
IM7/SP500 samples were subjected to the laser source for 
durations of 0.5 and 1.0 s. The IM7/PEEK sample was sub- 
jected to the laser source for 0.5 s. 

In this study. Lamb wave velocities were measured us- 
ing the same technique described in the second fatigue study. 
The velocity of the lowest order symmetric Lamb mode was 
measured in the 0° direction in both the undamaged and 
damaged regions of each of the samples. The velocities as 
well as the damage level for each of the different samples are 
shown in Table I. The listed velocities are the average of 
three measurements performed in each region. For each of 
the three measurements in a particular region, the values of 
the velocity varied less by less than 5% from the average. 

As can be seen from the data, the velocity actually in- 
creases in the area which was subjected to the laser source. 
This is contrary to what was seen in the thermally damaged 
samples discussed earlier in this paper. The main difference 
between the laser damage and the thermal damage previ- 
ously mentioned is that the laser subjects the sample to an 
extremely rapid thermal cycle. This type of heating could 


TABLE I. Data for laser damaged samples. 


Sample 

Irradiation 

(W/cm 2 ) 

Time 

(s) 

Fluence 

(J/cm 2 ) 

Undamaged 
Vel. (km/s) 

Damaged 
Vel. (km/s) 

IM7/SP500 

100 

0.5 

50.0 

6.01 

6.10 

IM7/SP500 

100 

1.0 

100.0 

5.91 

6.33 

IM7/PEEK 

100 

0.5 

50.0 

6.50 

6.56 


actually cause embrittlement of the sample, which would 
serve to increase the velocity of the Lamb wave in that re- 
gion. 

V. CONCLUSION 

The Lamb wave velocity measurements in this study 
were conducted at long wavelengths. This was done for sev- 
eral reasons. First, if the wavelength is large compared to the 
diameter of the fibers, composites can be treated as homoge- 
neous. Second, only the lowest order modes propagate in the 
frequency range where the contact measurements were con- 
ducted. Finally, the effective elastic properties of the mate- 
rial can be measured because local defects, which scatter 
high frequency waves, will not be observed at long wave- 
lengths. 

Lamb waves offer a useful technique for characterizing 
many forms of damage in composite materials. Lamb wave 
velocity measurements are better than conventional ultra- 
sonic measurement schemes (i.e., through-the-thickness 
measurements) because they can propagate over long dis- 
tances and are sensitive to the desired in-plane elastic prop- 
erties of the material. The propagation of Lamb waves de- 
pends on a variety of material properties: elastic stiffness 
constants, density, and thickness. As a material is damaged 
in some fashion, one or more of these material properties are 
altered. Since the Lamb wave velocity is directly related to 
these parameters, an effective tool exists to monitor damage 
in composites by measuring the velocity of these waves. 

The Lamb wave velocity is a quantitative measurement 
and it has been shown by this work to be an effective tool in 
monitoring both fatigue and thermal damage in composites. 
Thus the Lamb wave method can be used to verify the in- 
tegrity of a composite structure over the service life of a 
component. This is an important measurement for flight 
qualified composite materials that may have undergone some 
measure of nonvisible damage, which could seriously de- 
grade the performance of the structure. Since the Lamb wave 
velocity depends on a variety of material properties, an ideal 
technique exists to monitor composites as damage is in- 
curred. With the continued development of damage assess- 
ment techniques such as the Lamb wave method, the safety 
of such structures can be assured. 
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